Spinodal decomposition during the hadronization stage at RHIC? 
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The expansion of strongly interacting matter formed in high-energy nuclear collisions drives the 
system through the region of phase coexistence. The present study examines the associated spinodal 
instability and finds that the degree of amplification may be sufficient to raise the prospect of using 
the spinodal pattern formation as a diagnostic tool for probing the hadronization phase transition. 



PACS numbers: 25.75.-q, 64.70.-p, 47.75. +f, 47.54.- 



m 
o 
o 

(N 

<: 
c- : 



■ 

(N ■ 

oo : 
o . 

m ■ 

o ; 

^ : 

a: 
i 

D 

^ ■ 



X 



I. INTRODUCTION 

The central goal of high-energy nuclear collisions is to 
explore the expected phase transition from the familiar 
hadronic world to a plasma of quarks and gluons. This 
phase change of strongly interacting matter is expected 
to be of first order with a critical temperature T c in the 
range of 150 — 200 MeV, although this remains to deter- 
mined experimentally. 

When two heavy nuclei collide at sufficiently high en- 
ergies, such as those provided by the Relativistic Heavy 
Ion Collider (RHIC), the matter in the central rapid- 
ity region is formed well above the critical temperature 
T c in a state of rapid longitudinal expansion. The ex- 
pansion then causes the energy density e to continually 
drop, thus forcing the system through the critical phase 
region towards hadronization and subsequent chemical 
and kinetic freezeout. During this evolution the system 
maintains local equilibrium and one may employ fluid 
dynamics which has the advantage that the equation of 
state, p(e), enters explicitly in the equations of motion. 

The starting point for the present discussion is the ob- 
servation that the occurrence of a first-order phase tran- 
sition is intimately linked with phase coexistence and as- 
sociated spinodal instability. Indeed, a first-order phase 
transition occurs when the appropriate thermodynamic 
potential exhibits a convex anomaly . There is then an 
interval of energy density throughout which the pressure 
has a negative derivative, dp/de < 0. This anomalous be- 
havior identifies the region of spinodal instability where 
small deviations from uniformity are amplified, as the 
system seeks to break its uniformity and separate into 
the two coexisting phases consistent with the given s. 

However, in the absence of reliable models for the tran- 
sition region, fluid dynamical studies of nuclear collisions 
have usually employed the (constant) pressure describ- 
ing the phase-separated system, thus suppressingentirely 
the spinodal instability (see, for example, Refs. 0,SH) 
Employing a simple spline technique to obtain the ap- 
proximate form of the anomalous pressure curve, the 
present study examines the significance of this inherent 
instability and seeks to ascertain the degree to which it 
may generate spatial patterns that might offer a diagnos- 
tic tool for probing the phase transition. 



II. EQUATION OF STATE 

To estimate the behavior of p(e) through the region of 
phase coexistence, we employ a cubic spline function to 
interpolate between the hadron and plasma phases. As 
is commonly done 0, 0, 13 , we approximate the former 
phase by an ideal gas of hadrons, while we take a bag of 
non-interacting quarks and gluons for the latter. 

In general, the contributions to p and e from a non- 
degenerate particle of mass m are 
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where k = + for bosons and k = — for fermions. 

The hadronic phase is approximated as an ideal gas 
of 14 hadronic species i — n, K, rj, p, . . . , N, A, E, A with 
masses mi and degeneracies gi (which include antipar- 
ticles where appropriate). Since we are focussing on 
the central rapidity region of collisions at RHIC, we 
shall asssume that all the chemical potentials vanish, 
Ms = fJ-Q = fJ-s = 0- We then have 
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The hadronic equation of state is shown in Fig. ^ At the 
critical temperature, for which we use T c = 170 MeV, the 
pressure is p c = 54 MeV fnT 3 , while the energy density 
is £\ = 220 MeV fm -3 . The slope of p(e) gives the square 



of the sound speed at T c , {dp\ la d/de)\ 



0.27. 



For the plasma phase, we assume that the system can 
be approximated as an ideal gas of quarks and gluons 
confined within a bag of suitable (negative) pressure —B. 
The gluons as well as the u and d quarks are taken to be 
massless, m g , m q = 0, while the s quark is given a mass 
of m s = 150 MeV (though the results would change little 
if it were also zero) . The corresponding degeneracies are 



2 



= 16, g q = 24, and g s = 12. Then 

Pqgp(T) = [gg 
£q SP {T) = [gg 
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Thus the pressure starts at —B, while the energy den- 
sity starts at B. These quantities increase more rapidly 
than their hadronic counterparts and the plasma pres- 
sure overtakes the hadronic pressure at the critical tem- 
perature. This condition determines the bag constant, 
S 1 / 4 — 250 MeV. The resulting plasma equation of state 
Pqgp(s) is shown in Fig. The critical energy density of 
the plasma is e 2 = e qgp (T c ) = 2220 MeV fm~ 3 , so that 
the latent heat of the phase transition is Ae = e-i — e\ = 
2000 MeV fm~ . The function p qg p(e) is practically lin- 
ear (its slight non-linearity arises from m s ) and it has the 
slope (dp qgp /de) 2 —v\— 0.33 at T c . 
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FIG. 1: The equation of state for uniform strongly interacting 
matter (solid) as obtained by interpolating with a cubic spline 
through the coexistence region between an ideal hadronic gas 
for e < £i and a bag with quarks and gluons for e > E2- 
The dashed line through the phase coexistence region repre- 
sents the pressure of an equilibrium mixture of the two phases. 
The spinodal region extends between the two extrema of p(e) 
where uniform matter is mechanically unstable, dp/de < 0. 

Through the region of phase coexistence we employ a 
cubic spline function that matches the pressure and its 
derivative (i.e. the speed of sound) at both ends, 

p(e) = Pc + [fa - e)v\ - (e - Sl )v 2 2 ] l{e) , (6) 

where £ = (e— Sx)/Ae and£ = (e2~e)/ Ae. The resulting 
curve is shown in Fig. ^ This curve rises initially at E\ 
with the slope v\ but it gradually bends over and turns 
downwards. After crossing the mixed-phase line p = p c 
it starts bending upwards and finally joins the plasma 
curve at £2 with the slope v\ . The crossing occurs at the 
energy density e x given by (v\ + wf) £ x = v\e2 + v\s\. 

In the cubic spline approximation, the slope dp/de is 
parabolic in the phase coexistence region, 
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The inflection point in p(e) is at e = (si +e x +62)/% and 
its extrema (located at £a and eg) are separated by 



sb - £A = ^ [{£2 - ex) 2 + (e x - ex) 2 + (e 2 - £x) 
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The steepest slope of p(e) occurs at e and is given by 
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0.150 



( 9 ) 

This value is rather robust against changes in T c and m s . 
The £ dependence of the sound speed is shown in Fig. [21 
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FIG. 2: The magnitude of the sound speed, \v a \ = [dp/del 1 ^ 2 , 
as a function of e from zero to above the phase-coexistence 
region. While v s is approximately constant in each of the 
single-phase regions, it drops rapidly as the metastable phase- 
coexistence region is entered, and becomes zero at the spin- 
odal boundaries. It is imaginary inside the spinodal region 
where the thermodynamical potential is convex and sponta- 
neous phase separation occurs. In the mixed-phase scenario 
v 3 is zero throughout the entire phase coexistence region. 



III. ONSET OF SPINODAL DECOMPOSITION 

The system is regarded as an ideal fluid with no con- 
served charges. The energy-momentum tensor T^ v is 
then of ideal gas form , 



T" u = (e + p)u"u" - pg" 



(10) 



where = 7(1, v) is the four-velocity of the energy flow, 
and local conservation of four- momentum yields the basic 
equations of motion, d^T^ v = 0, i.e. 



d t [e + pYi 2 -p] = V • [(e + ph 2 v] 
d t [(e + p)i 2 v] = V[(e+p)-f 2 v 2 -f 
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These coupled equations are closed by the specification 
of the equation of state p(e). 

We now consider a small perturbations on a uniform 
static system, 



e(t,r) = eo + 5e(t,r) , Se(t, r) < e 



(13) 
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Since the flow speed v is of the order of Ss/eq, it may 
be regarded as small, v <C 1, and we may ignore terms 
of order v 2 and thus use 7=1. Insertion of e(t, r) from 
(JT3J into eqs. (|11I12|) then yields 

« d t E-V-[(e+p)v]&d t 5e-(e +po)V-v, (14) 
« d t [(e+p)v]- Vp^(e Q + P o)d t v-v 2 V5e, (15) 

where we have used Vp(e(r)) = (dpo/deo)^e = ufVe, 
with po = p( e o) being the pressure in the unperturbed 
uniform system. The speed of sound was displayed in 
Fig.0 Taking the gradient of (|15|) and using the derived 
relation l|14|) between 8e and v then yields the familiar 
equation for the propagation of sound waves, 
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For harmonic perturbations, 5e ~ exp(ifc-r — iuit), it 
yields a linear relation between the wave number and 
the frequency, u>k = v s k. Thus, if a physical system 
is brought into the spinodal region where dpo/deo < 
(such a quenching can be accomplished by rapid expan- 
sion and/or cooling), the frequency is imaginary and 
small deviations from uniformity evolve exponentially. 

The linear dispersion relation becomes unrealistic once 
the wave length of the perturbation, Afe = 27r/fc, has 
shrunk to the size of the range of the interaction be- 
tween the constituents of the fluid. At such small scales, 
the physical justification for employing fluid dynamics 
is questionable. However, a physically meaningful treat- 
ment can still be made within the framework of fluid 
dynamics by employing a suitable smoothing, as we shall 
now explain. 

The equation of state p(e) gives the pressure in a spa- 
tially uniform system as a function of its energy density. 
When the energy density varies sufficiently gently with 
position, one may obtain the corresponding local pres- 
sure as p(r) ps p(e(r)), i.e. one may calculate p as if the 
energy density had everywhere the same value as at r. 
This local density approximation becomes increasingly 
inaccurate as the scale of the spatial variations in e(r) is 
decreased, because the finite range of the forces reduces 
the response to fine ripples. A simple way to account 
approximately for the effect of the interaction range is to 
average the pressure over the interaction volume. While 
this procedure is in fact exact in certain cases (e^q. semi- 
classical one-body models of nuclear systems j(j), it is 
adopted here merely as a simple and transparent approx- 
imation which needs ultimately to be justified/improved 
on the basis on microscopic models of the system. 

Thus we assume that the local pressure p(r) can be ob- 
tained by a convolution with a (normalized) kernel g(r), 
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FIG. 3: The spinodal dispersion relation 119H at the energy 
density e where the pressure decreases most rapidly. Results 
are shown for three values of the range a. 



then proportional to the gradient of the energy density, 

Vp(r) « ^9kVe(r) = ikg k ^5e(r) , (18) 

where gk is the Fourier transform of g(r). Hence the re- 
finement amounts to simply modulating the local-density 
term dpo/deo in the dispersion relation by gk, yielding 
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As a result of this convenient factorized form, the largest 
collective frequency occurs at the same wave number ko 
for any energy density, thus enhancing the formation of 
a distinctive spinodal pattern of a characteristic scale. 

The spatial scale of the pattern is determined by the 
range of g(r), which one would expect to be in the range 
of 0.5 — 1.0 fm, with the precise value depending on the 
specific microscopic treatment. For simplicity, we employ 
here a gaussian form factor, g(r) ~ exp(— r 2 /2a 2 ), so we 
have gk = exp(— a 2 k 2 /2). The growth rates 7fc for the 
spinodal modes can then be written on a simple form, 
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This dispersion relation is illustrated in Fig. [3J 

As the system expands and the energy density tra- 
verses the spinodal region, the growth rate 7fc (£(*)) of a 
given mode increases from zero at £#, exhibits a max- 
imim at e, and then reverts to zero at ea- The to- 
tal growth factor is then approximately equal to Gk — 
exp(Ffc) 0, with the amplification coefficient Tk being 
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When the perturbation of the energy density is harmonic, 
Se(r) ~ exp(ifc ■ r), the gradient of the local pressure is 



The effective time during which spinodal instability is 
encountered is At e g « (tt/ A)(ts — £a) in the cubic spline 
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FIG. 4: The factor Gk = exp(Ffc) by which the amplitude 
of undulations of wave length Xk = 2n/k grows during the 
expansion through the spinodal region. 

approximation. For a quantitative estimate we use the 
value tB — tA ~ 4 fm/c gleaned from Ref. p|. The result- 
ing growth factors Gk are displayed in Fig. 0] 

After rapidly reaching a maximum, G(A) subsides rela- 
tively slowly, thus allowing some degree of amplification 
of longer wave lengths. If a = 0.5 fm, the maximum 
growth factor of G max w 8 is obtained for A « 2 fm, 
while a = 0.7 fm yields G max ss 4 for A « 3fm. The 
corresponding component of the density-density corre- 
lation function is proportional to the square of that, 
(fe(ri) exp(zri2 • k) fc(r 2 )) ~ G|, and can thus reach 
quite appreciable magnitudes. (Without the finite-range 
modulation factor gk there would be even larger amplifi- 
cation but no characteristic scale would emerge.) 

IV. DISCUSSION 

Spindoal decomposition is intimately linked to the oc- 
currence of a first-order phase transition and is important 
in many areas of physics as a mechansim for pattern for- 
mation. In nuclear physics it plays a key role for nuclear 
multifragmentation where it gives rise to highly non- 
statistical fragment size distributions with a preference 
for equal masses Q, H| ■ The present study is focussed on 
high-energy nuclear collisions and we have investigated 
the importance of spinodal decomposition while the ex- 
panding matter passes through the phase coexistence re- 
gion. Working within the framework of relativistic fluid 



dynamics and basing our analysis on commonly employed 
assumptions about the equation of state (supplemented 
with a simple spline procedure), we have found that a 
significant degree of amplification may occur while the 
matter passes through the mechanically unstable region 
of phase coexistence. 

Such an onset of spinodal decomposition may have in- 
teresting dynamical consequences, most notably the de- 
velopment of characteristic patterns in the energy den- 
sity, as the most rapidly amplified wavelengths grow pre- 
dominant. The spatial scale of any such emerging spin- 
odal pattern could in principle be probed by HBT in- 
terferometry. It is also conceivable that the azimuthal 
multipolarity of the favored modes would be reflected in 
the flow coefficients v m (of which vi quantifying ellip- 
tic flow is the best explored 0). Investigation of these 
prospects would clearly be interesting. 

If the appearance of spinodal patterns could in fact be 
examined experimentally, one could obtain fairly direct 
information on the equation of state in the critical region, 
since the relative prominence of various scales would re- 
flect the degree of spinodal instability encountered by the 
corresponding wave lengths. Conversely, if the absence 
of such pattens can be established it may be concluded 
that no spinodal region has been encountered and, con- 
sequently, the phase transition cannot be of first order. 

In view of the potential utility of spinodal decomposi- 
tion as a diagnostic tool, it would seem worthwhile to 
pursue this issue with more realistic dynamical treat- 
ments that take better account of the spatio-temporal 
features of the system (including both the rapid longi- 
tudinal expansion and the finite transverse extent which 
may affect the present simple estimates). In particular, 
the quantitative importance of spinodal decomposition 
could be elucidated with existing fluid dynamical models 
that have been modified to accomodate an equation of 
state containing the convexity anomaly which is a char- 
acteristic feature of first-order phase transitions. 
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